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C/3 . Abstract 

In a general two-scalar-doublet model without fermions, there is a unique source of 
1 _— I ■ CP violation, J\, in the gauge interactions of the scalars. It arises in the mixing of the 

three neutral physical scalars X\, X2 and A3. CP violation may be observed via different 
\ decay rates for X\ — > H + W~ and X\ — > H~W + (or, alternatively, for H + — ► X\W + and 

^NJ ■ ii -- — > X\W~ - depending on which decays are kinematically allowed). I compute the 

part of those CP-violating decay-rate differences which is proportional to J\. The CP- 
invariant final-state-interaction phase is provided by the absorptive parts of the one-loop 
OS ' diagrams. I check the gauge invariance of the whole calculation. 

a,: 
p<: 

CD . 

^ 1 Introduction 

> : 

There are general reasons for the interest in the possibility of CP violation in the scalar sector. 
CP violation is a necessary ingredient for the generation of the baryon asymmetry of the 
Universe |E|. It is believed that CP violation in the Kobayashi-Maskawa matrix is not large 
enough to explain that asymmetry 0. It has been speculated || that the scalar sector might 
provide the missing CP violation. 

There have been studies of possible signatures of CP violation in the scalar sector. It has 
been remarked ||] that the simultaneous presence of the three couplings Z°SiS2, Z°SiS 3 , and 
Z°S2S$, where Si, S2 and S3 are three neutral scalar fields in any model, implies CP violation. 
Similarly, the simultaneous presence of the three couplings SiZ°Z°, S2Z°Z°, and S1S2Z , 
represents CP violation. This is because the C quantum number of the Z° is — 1. Another work 
II has considered various CP-violating Lagrangians including scalars, fermions, and vector 
bosons, and has suggested looking for CP violation in the decay mode S — > Z°W + W~ 
occurring when, in the rest frame of the decaying neutral scalar S, the momentum distribution 
of the W + is not the same as the momentum distribution of the W~ - or, in a similar fashion, 
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in S — > Z°H + H~. The first of these CP- violating asymmetries has later been computed 
in the context of the two-Higgs-doublet model. However, the decay mode S — > Z°W + W~ is 
phase-space disfavored as compared to the simpler decay modes S — > W + W~ and S — > 
Other studies [0 have concentrated on CP-violating phenomena originating in the interplay of 
scalars and fermions, in particular the effects of top-quark physics. 

The aim of this work is the computation of a CP-violating asymmetry in the two-Higgs- 
doublet model without any fermions. The model has gauge symmetry SU(2)<S>U(1), which 
is spontaneously broken to the U(l) of electromagnetism by the vacuum expectation values 
(VEV's) of the two Higgs doublets. I look for CP violation involving solely the gauge interac- 
tions of the scalars. For simplicity, I do not consider the presence of fermions, which presence 
would lead to extra sources of CP violation, both in the fermion sector, and in the interplay of 
the fermion and the scalar sectors. I also omit possible sources of CP violation in the cubic and 
quartic interactions of the physical scalars. Those scalars are two charged particles H ± , with 
mass tuh, and three neutral particles X\, X 2 and X$, with masses mi, m 2 and 7713, respectively. 
Besides, the spectrum of the model includes the massive intermediate vector bosons W ± and 
Z°, with masses mw = 80 GeV and mz = 91 GeV respectively, and the massless photon. 
For a fairly large range of the masses of the scalars, either the two decays X\ — > H + W~ and 
X\ — » H~W + , or the two decays H + — > XxH 74 " and H — > XiW~, are kinematically allowed 
(the neutral scalars may be numbered so that X\ is the scalar for which one of these couples 
of decays is allowed). Then, the possibility of a CP- violating difference between the rate of 
one decay and the rate of its CP-conjugated decay exists. It is my purpose to calculate that 
difference. 

It has recently been observed || that the two-Higgs-doublet model has one and only one 
source of CP violation in the gauge interactions of the scalars. I describe it briefly. Because 
the U(l) of electromagnetism is preserved in the symmetry breaking, we can, without loss of 
generality, choose a basis for the two scalar doublets in which only one of them, Hi, has a 
VEV v, while the second one, H 2 , does not have a VEV. The two doublets in that basis can be 
written 

Hl= {v + iH + iG )/^) ' H2= ( (Ji + iOA/O ' (1) 

G + and G° are the Goldstone bosons, which become the longitudinal components of the W + 
and Z°, respectively. H, R, and I are linear combinations of the three neutral scalar fields 
Xi, X2, and X3, which are the eigenstates of mass. Those linear combinations are given by an 
orthogonal matrix T, 

( H\ ( Xi\ 

(2) 

Without loss of generality, we can assume T to have determinant +1. Then, the following 
useful identities follow: 

3 

TiiFzj — T^jT-a = ^ €ijkTik , (3) 
fc=i 

where is the totally antisymmetric tensor with €123 = +1. There is CP violation in the 
gauge interactions of the scalars || if and only if m 1; m 2 and m 3 are all different, and 

Ji = T n T 12 T 13 (4) 
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is non-zero. The quantity J\ has in the two-Higgs-doublet model a role analogous to the one of 
Jarlskog's |J J in the three-generation standard model. Notice however that, here, there are 
in principle other sources of CP violation, in the cubic and quartic interactions of the scalars. 

1 will neglect those extra sources of CP violation throughout this work. 

It is important to remark that, though J\ represents CP violation in the mixing of the 
three neutral scalars, this source of CP violation has nothing to do with the fermions and 
with the identification of, say, H and R as being scalars, and / as being a pseudoscalar. That 
identification can only be done when a specific Yukawa Lagrangian, coupling the two scalar 
doublets to the fermion sector, is considered, which I do not do here. Specifically, as is clear 
from Eq. |4], R and I play a completely equivalent role in J\ — indeed, as long as there are no 
Yukawa couplings, R and / may rotate into each other by a simple U(l) rephasing of H 2 . Also, 
Ji cannot be the source of, say, CP violation in the kaon system. If fermions are introduced in 
the model, the mixing of the neutral scalars will in principle lead to more CP violation than 
simply Ji, because of the Yukawa interactions of the scalars with the fermions. 

2 General features of the calculation 

Let us consider how CP violation proportional to J\ arises in the decay modes that I consider 
here. All the vertices needed for the calculations in this paper have been listed in Figure 1. 
The tree-level diagram with incoming particles W~, H + and X% is proportional to i(T 2 \ — iT 3 \), 
while the diagram with incoming particles W + , H~ and X\ is proportional to — i(T 2 i + iT 3 i). 
Now take a look at Figure 2, in which all the one-loop diagrams which lead to CP violation when 
interfering with the tree-level diagram are collected. Consider for instance the first diagram, 
with a loop of W + W~ , and then X 2 , as an intermediate state. That diagram is equal to 
2 6 T 11 T 12 (T 2 2 — iT 32 ), times a certain momentum integral The seven factors of % come, three 
from the vertices, three from the propagators, and one from the Wick rotation in the momentum 
integral. Therefore, the interference of this diagram with the tree-level one is proportional to 
the real part of (— i)(T 21 + iT 31 )i 6 T u T 12 (T 22 - iT 32 )iI k = T u T 12 (T u Ti 2 + iT 13 )I k . The T-matrix 
factor has an imaginary part equal to J\. Therefore, if the momentum integral has an absorptive 
{i.e., imaginary) part, then the interference term will include J\ times that absorptive part. 
This is CP-violating. The absorptive part of the integral plays in the calculation the role of a 
CP-invariant final-state-interaction phase, which allows J\ to manifest itself. 

We can check in a similar fashion that the absorptive parts of all other nine one-loop 
diagrams in Figure 2 lead, when one considers the interference of those diagrams with the 
tree-level one, to CP violation. Indeed, a careful study of the model and all its vertices shows 
that the ten diagrams in Figure 2 are the only ones which lead to CP violation proportional 
to J\ in this process.^ The CP violation manifests itself in a difference of the decay rates of 
X\ -> W+H- and X x -> W~H + , or of H + -> X{W + and H -> X X W~ , whichever pair of 
decays is kinematically allowed. 

Let me be more explicit. At tree- level, the amplitude for the decay X± — > W + H~, or for the 
decay H + — ► XiW + , is (e u P^)ig(T 2 i — iT 31 ), from the tree-level vertex. Here, e u is the polariza- 

1 There are other diagrams which may also lead to CP violation in this process, but which include other 
sources of CP violation, in the cubic scalar interactions. I neglect those extra sources of CP violation, just as I 
neglect fermionic sources of CP violation. 
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tion vector of the outgoing W + , and P^ is the incoming momentum of the H ± . This is because 
the polarization vector of the W + is orthogonal to the momentum of that vector boson. At 
one-loop level, for the same reason, each diagram in Figure 2 contributes M*. = (e u P^)g 3 CkiIk- 
Here, C& are the various % factors and T-matrix factors from the vertices and propagators in 
the diagram, and ilk is the momentum integral, with the % coming from the Wick rotation. 
The amplitudes for the CP-conjugated decays are, at tree-level, {e ll Pj I ){—i)g{T2i + iT^i), and, 
at one-loop level from Figure 2, (e u P^)g 3 (—Cl)Hk- Notice that, while the momentum integral 
is the same, the vertex factors are complex-conjugated. Then, the CP- violating asymmetries 
are 



or 



BR(Xx -> W + H-) - BR(X 1 -> W~H+) 
BR(Xi -> W+H-) + BR(Xi W-H+) 

BR(H+ -> X 1 W + ) - BR(H~ -> X X W~ 



BR(H+ -> X 1 W+) + BR(H~ -> XiW" 



2 1° Im[(T 21 +^T 31 )C fc ]Re(zJ fcj _ 

^ 1 ^2 ' W 



fc=l 



1-2?! 



where I used the orthogonality of T to write T 2 2 x + T 3 2 x = 1 — T n . The imaginary part of 
(T21 + iTzi)Ck is simply Ji times a number — typically ±1 or ±2. The momentum integral 
ilk (the i is from the Wick rotation) has a real part if cuts in the corresponding diagram lead 
to absorptive parts. Notice that in Eq. [5] I have used the approximation of taking, in the 
denominator, only the square of the modulus of the tree-level contribution to the amplitude. 
It is clear from Eq. [5] that the asymmetry will be of the form 

2 2T11T12T13 

9 — j W~ A ' ( 6 ) 

1 1 11 

where A represents the sum of the absorptive parts of all the diagrams in Figure 2, weighted 
by appropriate numbers ±1 or ±2 (see the preceding paragraph). As one is interested in how 
large the asymmetry can be, I now consider the mixing-matrix factor in Eq. |^. Because one is 
constrained by the orthogonality condition 2\ 2 2 + T 2 3 = 1 — T^, it is obvious that that factor 
will be maximal when |2i2| = | ^Z~i3 1 , and we will then have 

22iiTi 2 2i3 _ . . 

l-Tl ~ Tu ■ (7) 

Clearly, the order of magnitude of this quantity Tu is 1. It is at this point important to 
remark that, in this specific example, the CP-violating asymmetry approaches its maximum 
when the decay rate decreases. Indeed, as T n — > 1, the asymmetry becomes potentially larger 
(as long as |T 2 i| remains equal to |T 31 |), but the decay rate, which is proportional to 1 — T^, 
approaches zero. Similarly, the decay rate becomes larger if Tn — > 0, but then J x — > and 
the CP asymmetry also vanishes. This situation is reminiscent of the case of CP-violating 
asymmetries in decay modes of the B° mesons, which are generally predicted to be larger when 
the branching ratios are smaller, and vice-versa. 



3 Gauge invariance 

A way to check that the diagrams in Figure 2 are all the relevant ones is to check whether their 
computation yields a gauge-invariant result. Because I just want to compute the absorptive 
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parts of the diagrams, which are finite, it is sufficient to compute the diagrams in the unitary 
gauge, in which no Goldstone bosons and no ghosts are present. However, in a general 't-Hooft 
gauge, in each of the diagrams in Figure 2, a G ± can be used instead of the W ± , or a G° 
can be used instead of the Z°. In the two diagrams which have a loop only of W ± or of Z Q 
(diagrams 1 and 2), ghost loops must also be considered in an 't-Hooft gauge. Now, in an 
't-Hooft gauge, the W ± propagator contains an extra piece (relative to the unitary gauge) in 
which the W ± has an unphysical squared mass W . Similarly, the charged Goldstone bosons G ± 
and the charged ghosts c ± have unphysical squared mass W in an 't-Hooft gauge. When W is 
not infinite as in the unitary gauge, each diagram by itself contains a H^-dependent absorptive 
part. However, all those unphysical absorptive parts must cancel out when one considers the 
whole set of diagrams. The same thing can be said about the propagators of the Z°, which has 
a piece with unphysical squared mass Z, and of the Goldstone boson G° and ghost c°, which 
have squared mass Z. (In principle, Z ^ W '.) The whole set of propagators is given in Figure 
3. The sum over all diagrams of all the absorptive parts must be independent of both W and 
Z. I have checked that independence. 

To be sure, gauge-independence only applies to an observable quantity. Thus, gauge- 
independence in this case only occurs when 1) the three external particles are all on mass 
shell, that is, P\ = m 2 H , P^ = m 2 w , and P 2 — in 2 , where Pjj, Pyj and P\ are the incoming 
momenta of the external H ± , W ± , and X ± , respectively; 2) one suppresses from the amplitude 
all terms proportional to Pfo, because the amplitude must be multiplied by the polarization 
vector e u of the W ± , and e v P^ = 0; 3) one considers, from each one-loop diagram, only the 
part which is proportional to J\ upon interference with the tree-level amplitude; 4) one only 
considers the absorptive part of each one-loop diagram. One has also to take into consideration 
that the gauge- dependent absorptive parts sometimes cancel between two similar diagrams with 
intermediate virtual particle X 3 instead of X 2 , whenever those absorptive parts do not depend 
on the mass of that intermediate particle (m 3 or m 2 ); while other absorptive parts cancel among 
different-looking diagrams. 

4 Contribution of each diagram 

I present in this section the results of the computation of each of the ten diagrams in Figure 2. 

I am only interested in the absorptive part of each diagram. Those absorptive parts are 
simpler to compute when the power of the integration momentum k in the denominator is 
lower; in particular, if only one propagator k 2 — m 2 occurs in the denominator, there is no 
absorptive part. Moreover, I want to display explicitly how the cancellations which lead to 
gauge invariance of the whole result occur. Therefore, my method has been the following. First 
I have added, for each particular diagram, the contributions with gauge bosons, the ones with 
Goldstone bosons, and the ones with ghosts. Then, I have simplified as much as possible the 
numerators of the integrands, in such a way that factors in some terms of those numerators 
cancel out some of the propagators in the denominators, thereby reducing the overall power of 
the integration momentum k in the denominator. As a result, each diagram becomes the sum 
of three parts: one piece with an overall k 2 in the denominator, leading to no absorptive part; 
one piece with k 4 in the denominator; and — for the diagrams 7 through 10 only — one piece 
with k 6 in the denominator. I write down here the two latter pieces only, because they are the 
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only ones which lead to absorptive parts in the integrals. 

For each diagram, I have taken separately the factors with T-matrix elements, and all the 
i factors from both the vertices and the propagators, multiplied that by the factor (T 21 + iT 3 i) 
arising in the interference with the tree- level diagram, and taken the imaginary part of the result. 
That imaginary part is always a multiple of J\. This I did because it is only the quantities 
Im[(T2i + iT^i)Ck\ which are relevant for the computation of the asymmetry, according to 
Eq. [|. Also, it is only for these quantities Im[(T 2 i +iT 31 )Cfc] that the cancellation of the gauge 
dependence must occur. As a consequence, each momentum integral presented below should be 
looked upon as being an ilk in the notation of section 2 (the i arising when the Wick rotation 
is performed), and only the absorptive parts of each such momentum integral are meaningful 
in this context. 

Each diagram in Figure 2 is in reality two diagrams, one in which an intermediate state 
A 2 appears, and another one in which the intermediate state is A3. Those two diagrams have, 
as is easily seen, opposite signs of J\ in the interference term Im[(T 2 i + iT 3 i)Ck]- Therefore, 
if the corresponding momentum integrals happen not to depend on the masses m 2 or 777.3 ~ 
because the propagator of A 2 or of A3 might have been cancelled out by a similar factor in the 
numerator of the integrand — then those terms become irrelevant for the present computation. 
As a consequence, for each diagram, I only write down those integrals which do depend on m 2 
or m 3 . 

As before, Pi, Pw and Ph are the incoming momenta of the scalar X\, the gauge boson W~ , 
and the charged scalar H + . I have systematically set all three external particles on mass-shell, 
and used whenever possible e u P w = in order to simplify the amplitudes. 

I first define the various combinations of the integration momentum k and of the incoming 
momenta, and of the masses, which appear in the denominators of the momentum integrals. 
They are 



D 



D 2 



k 2 — m\ 
k 2 + 2k 
k 2 + 2k 
k 2 -2k 
k 2 + 2k 
k 2 + 2k 




P w - 2k ■ P H - m\ + 2m 2 w + m\ 



(8) 
(9) 
(10) 




(12) 
(13) 



Notice that only D 4 depends on m 2 . 

It is convenient to define the "triangular function" A: 



X(A, B, C) = A 2 + B 2 + C 2 - 2(AB + AC + BC) . 



(14) 



This function is negative if and only if one can form a triangle with sides of length \f~A, \^B, 
and \/C. 

I now present the results for the first three diagrams. Diagram 1: 



Diagram 2: 

M 2 



m% 



rrin 



x 



e P u n 3 T - 

(mf — m 2 )(mf 
—m\ + Am\m 2 z — 12m| 
8m| 



-ml) 

ci 4 fc 1 

(27r) 4 Di(m w — > m z )D 2 (m w — > m z ) ' 



(16) 



(There is an extra factor 1/2 in diagram 2 relative to diagram 1, which is a symmetry factor, 
due to the two identical Z°'s in the loop.) Diagram 3: 



M, 



m 



\(ml,ml v ,m 2 H ) f d 4 k 



(ml 



ml) (ml — m 2 ) 



2m 2 w 



r crt 
J (27 



(2tt) 4 DiD 3 ' 



(17) 



The results for each of these three diagrams are separately gauge-invariant. This is because the 
gauge-dependent piece in each of them (the piece depending on either W or Z) is the same, for 
each case, in the diagram with an intermediate X 2 , and in the diagram with an intermediate 
X 3 . As the sign of the J\ factor is opposite, the gauge-dependent piece cancels out between 
the diagrams with intermediate X 2 and those with intermediate X 3 . The same phenomenon 
partially occurs in all other diagrams but, in each of them individually, some gauge dependence 
always remains, as is seen in the following. 
Diagram 4: 



M A 



d A k 



^(m 2 H 



m 2 2 ) 



(2nY 1 
\(m\,m 2 w ,m\) 



+ 



k v {mw + m\ 



m 2 2 ) 



H AD x D 4 m 2 w {m 2 H - m 2 w ) 



4m 2 w D 4 (ki - W) 2D 1 D 4 (m 2 H - m 2 w ) 
m 3 ) \ . 



Am\y 



\m 2 



(18) 



Diagram 5: 



M R = 



+ 




m z - ?7% 



m 



w 



J (2-) 4 r 



2 2 

m 2 — m\ 



ml + m\ 



Diagram 6: 



M 6 = t v P v H g\h X - J 



+ 



d 4 k 

\(m 2 z ,m 2 ,ml) 
DAw^K - ml) 



4m| 



m% — m 



Am 2 z D 4 [(k + P w ) 2 - Z] 
m 3 ) \ . 



im 2 



(19) 



m 



w 



D 4 [(k 
- (m 2 <-> m 3 ) \ . 



Pw) 2 ~ Z\ 



(20) 



Diagram 7: 

M 7 = 

+ 



2m w 



m 



Am 2 w 



1 2 2 2 

i m 2 — m\ — m z 



+ 



D 4 D 6 
(m 2 <-> m 3 )} . 



m 2 z 



d A k 

1 

DaD* 



(Ph - kf 



(m\ — ml) 



m 2 z D 4 [(k + P w ) 2 - Z] 
+ (2m w + m| — 2m 2 H — m\ — m 2 2 ) 



1 



D 4 D 5 D 6 



(21) 
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Diagram 8: 



e v g 3 Ji 



2m 2 w 



d A k 



{k - PhY 



H 



+ 



m 



H 



m 



W 



D X D A 



+ 



w 



+ 



1 



D 3 D 4 D l D 3 D A 



D 4 (k 2 - W) 

to 4 4,2 2 2 2 

{om w — m H + m l m H — m 1 m 2 



,22 22 22 o 2 

+m H m 2 — m x m w — m 2 m w — 2m v 



w m H 



(m 2 <-> m 3 ) 



(22) 



Diagram 9: 

M 9 = 
+ 
+ 
+ 

Diagram 10: 



-tv g 3 Ji 



d A k \ (k - PhY 



[m 



H 



{k-P u 



Hj 



4 
1 



(2vr) 4 
1 



+ 



{k-P H y {m 2 H -ml) 



Am 2 w D A [(k + P x ) 2 - W\ Am 2 w D A (k 2 - W) 



(k 



PhY m 2 



n 



D 2 D 4 D ± D 4 
Qmly + 2m 2 H + m\ 



D 1 D 2 D 4 
m\(m 2 H — ml) 



pv 



Am 2 w 



mw 
2m 2 . 



DoDa 



+ 



DiD 4 



m 



1 UV 



Am 2 w 



(k - PhY 



- (m 2 



m 3j 



(23) 



Mi 



10 



3 1 /■ C? 4 fe f 1 

~ 6 " 9 1 4 J (2vr) 4 { (fc 2 - W) [(k + P w ) 2 



m\ 



Z] 



in 



D A [(k + P w ] 



- k v 
+ k v 



+ 



2 -Z] 
1 



(k - 2P H 



+ k" 



mr — m% 



m% 



D 4 (k 2 - W) 



{k - 2P H Y 



m 



D X D 4 D 4 DJ D ± D 4 D 5 
(m\ — m 2 H )(m\ — m%) 1 

m 2 w P> X D 4 D*> 



, ry ry rt rt 

(4m 2 — 2m w + 2m l + 2m H — m 



k v 



mr, 



mlyDxD^ 



+ {2P H -kY 



3m 2 , + ^ m w — ml — 2m 2 H 



D 1 D 4 D 5 



+ 



mi — m 2 



m% 



D A D 5 



+ 



m 2 H — ml 



1 



m\ 



m 



w 



D 1 D 4 DiD A D z 



(m 2 -> m 3 ) 



(24) 



It is simple to check now that most terms depending on either of the unphysical squared- 
masses W or Z cancel among M4 through M\q. The only two exceptions are the first terms in 
the curly brackets in the expressions for Mg and for M\q. The term in Mg contains the integral 
of (k - P H ) U / D 4 [(k + P^ 2 - W], while the one in M w contains the integral of k u / (k 2 - W)[(k + 
Pw) 2 — Z]. However, it is easily found that the absorptive parts of both these integrals are 
proportional to P-^, and therefore give a vanishing contribution to the amplitude. This ensures 
the gauge invariance of the whole computation. In each of Eqs. [L8] to [24], therefore, one only 
has to consider only the terms independent of W and of Z. Those terms are the ones that 
would have been obtained had the computation been performed directly in the unitary gauge. 
Some of those terms, however, still yield zero absorptive contributions. 
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5 Absorptive parts 



In this section I briefly review the computation of the absorptive parts of the integrals. 
Consider first an integral with k 4 in the denominator: 

4 J (2tt) 4 \{k + P X Y - m\\ [(k + P y ) 2 - m|] ' 1 ° J 

where r and s are some coefficients. Let (Px — Py) 2 = m c > 0- The integral is divergent. 
After introducing a Feynman parameter x, and using dimensional regularization, one obtains 



ih = Jo dx [r(Py ~ Px)Vx ~ rP v + sP ^ ( 



-7-ln— ~ 2 , (26) 



4 - d 4vr/i 



where 

A = m c x 2 + (m^ — m^ — m 2 c )x + m^ . (27) 

The absorptive part arises when A is negative for part of the integration domain of x. This 
happens when mc is larger than + mg. Substituting ln(A/47r/i 2 ) by — in and integrating 
over the part of the interval [0, 1] in which A is negative, one gets 

D ,._ , 1 v /A ( m ^' m B' m c) 

Re ui 4 ) = 



x 



167r m c 



m 



+ *P£ . (28) 



Using this result, one sees that the absorptive parts of some of the integrals in the previous 
section are proportional to P w , and therefore irrelevant. This is true not only of the two 
gauge-dependent integrals which did not cancel out, as I pointed out in the last paragraph 
of the previous section; it is also true of the absorptive parts of the integrals of k l/ /(D 1 D 5 ) 
(which arises in M 10 ), of (k - P H ) U /(D 2 D 4 ) (in M 9 ), of (k - P H ) U /(D 3 D 4 ) (in M 8 ), and of 
(k — Ph) u I {DaDq) (in M-?). Eliminating all of these, it is seen that the only remaining integrals 
of the type J 4 either have denominator D\D\, or have denominator P4P5. The first ones have 
an absorptive part if m# > mw + m 2 , and therefore H ± can decay to W ± X 2 , and the second 
ones have an absorptive part if m x > m z + m 2 , i.e., if X 1 — > Z°X 2 is possible. 

Now consider the integrals with k 6 in the denominator. In our problem, they are all of the 
form 

r d 4 k e v (rW + sP&) 

"6 = / 77^17170 -2TT7I , 137771 TTT^ I3T > ( 29 ) 



(2tt) 4 (P - mi) [(£; + Pi) 2 - ml] [(k - P H ] 



m\ 



c\ 



where r and s are coefficients. Specifically, for the integral of 1/(P 4 P 5 P 6 ) present in diagram 7 
we have = mz, = m 2 and mc = m#; for the integral of 1/ (D1D3D4) present in diagram 
8 we have m^ = mw, m# = m# and mc = m 2 ; for the integral of \j{D\D 2 D±) present in 
diagram 9 it is m^ = m# = mw, and mc = m 2 , and, for the integral of 1 / (D1D4D5) present in 
diagram 10, it is m^ = m 2 , m# = mz and mc = miy. The integral Iq is finite and, eliminating 
its part proportional to P w , one obtains 



9 



where 



A = A Y y 2 + B Y y + m\ , (31) 
A Y = m^x 2 + (m? - m 2 H - m\)x + m 2 H , (32) 
By = (m^ — m 2 + m^ — m 2 c )x + m^ — m 2 A — m 2 H , (33) 

and e + in a positive infinitesimal quantity. Using now the fact that the imaginary part of 
1/(A — ie + ) is z7r5(A), and then integrating over the Feynman parameter y, thereby getting rid 
of the Dirac S function, one gets 

B Y m\\ 1 



Re(z/ 6 ) = ^ dx [s— -+r-^-2r-^\ . . (34) 

Here, the integral must be performed over that part of the interval [0, 1] for which the two 
values of y defined by the condition A = are positive and smaller than 1. That sub- interval 
of [0, 1] must be found for each of the four integrals of the type Iq individually It is larger or 
smaller according to whether some cuts can be made in the corresponding Feynman diagram: 
namely, whether m# > rriA + mc or not, and whether mi > ttia + me or not. In any case, the 
integral over x can be performed analytically, and it yields a rather complicated logarithmic 
function. 



6 Results 

We now have all the ingredients needed to compute the asymmetry. The asymmetry is equal 
to g 2 ~ 0.43 times a mixing-matrix factor, studied in section 2, which should be of order 0.1 to 
1, times the sum A of all absorptive parts, which itself includes a suppression factor l/(167r). 

Now, one should note that the absorptive parts of diagrams 1, 2, 3 and 6 all diverge when 
m2 (or tus) approach mi. This is simply because in those diagrams one has a scalar X 2 (or 
X 3 ) propagating with momentum Pi such that P 2 = m\. Those divergences do not cancel in 
the absorptive parts, because the specific values of those absorptive parts depend on different 
parameters: in diagram 1, on the decay width X\ W + W~, in diagram 2, on the decay 
width Xi — > Z°Z°, in diagram 3, on the decay width X\ — > H + W~, and, in diagram 6, on 
the decay width Xi — > Z°X 2 or Z°X 3 . Of course, we know that these divergences are not 
genuine, they might be eliminated by a proper treatment in which one would take into account 
the finite width of the propagating X2 or X3, and besides, from a different line of reasoning 
||, one knows anyway that CP violation disappears and J± loses its meaning when the masses 
of any two of the three neutral scalars become equal. A proper treatment of these divergences 
at mi = m 2 or mi = m.3 would lead me far astray, and therefore I simply avoided, in general, 
considering the region of the parameter space in which either m 2 or are close to mi. 

Avoiding those regions in which the present approximation loses its validity, I find that the 
sum A of all absorptive parts is typically of order of magnitude 10~ 3 or 10~ 2 . A few examples 
are presented in the following table. 
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i is worth remarking that the total absorptive part A is always the final result of substan- 
tial cancellations among the absorptive parts, with different signs, of the various individual 
diagrams. 



7 Conclusions 

In this paper I have presented a model calculation of a CP-violating asymmetry in the two- 
scalar-doublet model. The asymmetry chosen has been the different decay rates for Xi — > 
H + W~ and for X x — > H~W + . I believe this to be a quite interesting place to look for CP 
violation in the scalar sector, even if the present calculation turns out not to be very relevant. 
This might happen mainly because I have taken into account only one source of CP violation, 
Ji in the gauge interactions of the scalars, while I neglected further sources of CP violation, in 
the cubic scalar interactions and in the Yukawa interactions with the fermions. Even within the 
context of the approximation that I have used, in which the decay proceeds essentially because 
of a tree-level vertex, while the CP violation arises because of the interference of that tree-level 
vertex with one-loop diagrams with absorptive parts, those further sources of CP violation will 
in principle lead to extra contributions to the total asymmetry. 

My interest here has been to ilustrate the specific way in which J\ arises in the computation 
of a CP asymmetry. I observed that there is a kind of balance between the CP asymmetry and 
the decay rate in this specific case — but only if the only source of CP violation is taken to be 
Ji -, with a large asymmetry being possible only when the decay rate is small, and vice-versa. 

I found that the asymmetry can attain values of order 10~ 2 . These values would increase or 
decrease if the interference with other sources of CP violation in this mode were constructive 
or destructive. 

Because of the presence of gauge bosons in the internal lines of the one-loop diagrams that 
I had to compute, I found it convenient to check the gauge invariance of the whole calculation. 
I checked that the fictitious masses that appear in the propagators of the W ± and of the Z°, 
and of the corresponding Goldstone bosons and ghosts, in a general 't Hooft gauge, lead to 
gauge-dependent absorptive parts for the individual diagrams, which however cancel out when 
all the diagrams which lead to CP violation proportional to Ji are considered. This constitutes 
a good check that one did not omit any diagram. 

I thank Lincoln Wolfenstein for suggesting to me the computation in this work, and for 
reading the manuscript. Ling-Fong Li and Martin Savage made useful technical suggestions. 
This work was initiated during visits to the Stanford Linear Accelerator Center, and to the 
Institute for Nuclear Theory at the University of Washington, the hospitality of which I ac- 
knowledge. This work was supported by the United States Department of Energy, under the 
contract DE-FG02-91ER-40682. 
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FIGURE CAPTIONS 

Figure 1: The vertices needed for the computation of the diagrams in Figure 2, and the 
Feynman rules for their values. All the particles are incoming particles, and all the momenta 
are incoming momenta. The indices k, j and / may be 1, 2 or 3. The meaning of the parameters 
W and Z is explained in the caption of figure 3. Other vertices which might at first sight be 
relevant, like W~W + G° and H~H + G°, do not exist. 

Figure 2: The one- loop diagrams with external incoming particles W~ , H + and X±, which 
lead to CP violation upon interfering with the tree-level diagram, if the integrals have absorptive 
parts. In each case, the W ± and the Z° in internal lines may be substituted by the corresponding 
Goldstone bosons G ± and G°, respectively; in the first two diagrams, they may be substituted 
by ghost loops as well. 

Figure 3: The propagators relevant for the computation of the diagrams in Figure 2. W 
and Z are unphysical parameters with dimension of squared mass, which also arise in the ghost 
vertices in Figure 1. The final physical results are independent of both W and Z. 
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